Crossover from skin-mode to proton-neutron-mode in El excitations of neutron-rich 

nuclei 



m 
o 

(N 



o 

=5 

(N 
> 

in 
o 
o 



(N 



X 



H. Nakada/'B T. Inakura, 1 * 2 and H. Sawai 1 

1 Department of Physics, Graduate School of Science, Chiba University 
Yayoi-cho 1-33, Inage, Chiba 263-8522, Japan 
2 Nishina Center for Accelerator-Based Science, RIKEN 
Hirosawa 2-1, Wako, Saitama 351-0198, Japan 
(Dated: February 7, 2013) 

Character of the low-energy El excitations is investigated by analyzing transition densities ob- 
tained from the RPA calculations in the doubly-magic nuclei. We propose a decomposition method 
of the El excitations into the pn mode [i.e. oscillation between protons and neutrons) and the skin 
mode (i.e. oscillation of the neutron skin against the inner core) via the transition densities, by 
which their mixing is handled in a straightforward manner. Crossover behavior of the El excita- 
tions is found, from the skin mode at low energy to the pn mode at higher energy. The ratio of 
the skin-mode strength to the full strength turns out to be insensitive to nuclide and to effective 
interactions in the energy region of the crossover. Depending on the excitation energy, the observed 
low-energy El excitations are not necessarily dominated by the skin mode, as exemplified for 90 Zr. 
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Introduction. While it has been known that the El 
excitations in atomic nuclei are predominantly concen- 
trated in the giant dipole resonance (GDR), recent ex- 
periments have disclosed that sizable El strengths are 
present at lower energies in a number of TV > Z nu- 
clei; in nuclei far off the j3 stability such as 68 Ni [l[ 
and 130 ' 132 Sn Q, and in medium- and heavy-mass nu- 
clei along the /3 stability line, 90 Zr J3L JjK the Pb iso- 
topes [5j-|9| and the N = 82 isotones [1CH17| . These low- 
energy El strengths are often called pygmy dipole reso- 
nance (PDR) , though their physical content yet remains 
obscure. The GDR is a typical collective excitation of 
nuclei, in which the constituent protons oscillate against 
neutrons {pn mode), keeping the center-of-mass (cm.) 
at rest. A possible interpretation of the PDR is an os- 
cillation of the neutron skin against the inner core (skin 
mode), as su gges ted in pioneering studies with macro- 
scopic models [li-iO]. However, recent calculations using 
the random-phase approximation (RPA) do and do not 
support the skin-mode picture of the observed PDRs, de- 
pending on nuclides as well as on effective interactions or 
energy density functionals. A systematic calculation 21 1 
shows that the low-energy El strengths well correlate to 
the neutron skin thickness in some regions of nuclides, 
but not in other regions. The low-energy El strengths 
are quite relevant to reaction rates in nucleosynthesis. 
Moreover, the low-energy El strengths could correlate 
to density-dependence of the nuclear symmetry energy, 
and therefore to structure of neutron stars. It is impor- 
tant to comprehend energy- and nucleus-dependence of 
the low-energy El excitations, to which inspection of its 
character should be crucial. 

The El transitions have isovector (IV) character in the 
long wavelength limit. In the pure pn mode the transi- 
tion densities should be IV at any position. However, the 



transition densities at each position may have isoscalar 
(IS) character in the skin mode, although the IS com- 
ponent should vanish after integration (see Eq. (|5|) given 
later). In Refs. [l3|, 15, it has been shown experi- 
mentally that low-energy 1~ states in 138 Ba, 140 Ce and 
124 Sn can be excited by (a, a'7) as well as by (7,7'), m 
contrast to higher-energy states excited only via (7, 7'). 
This seems to support the skin-mode interpretation of 
the low-energy El strengths in these nuclei. The isospin 
character of the low-energy transitions has been investi- 



gated for 140 Ce [22 1 



68 Ni, 132 Sn and 208 Pb [23, |2j by 
the RPA calculations. It has been pointed out in those 
studies that character of the transitions varies as energy 
increases. Similar results have been reported from calcu- 
lations in the quasiparticle-phonon model [25j. However, 
in these experimental and theoretical studies, character 
of the El transitions does not necessarily switch from the 
skin mode to the pn mode sharply. It is hence desired 
to investigate character of the low-energy El transitions 
with taking admixture of the two modes into consider- 
ation. Of particular interest is how the degree of the 
mixing varies for increasing energy and/or mass number. 
Since the transition operators to l - states have different 
radial-dependence between the IS and the IV channels 
so as not to contain the cm. motion, the IS and IV 
transition strengths cannot directly be compared. Wo 
should also notice that the oscillation of the neutron skin 
against the inner core should provide transition density 
dominated by neutrons outside the core, which has IV as 
much as IS characters. It is not easy to draw quantita- 
tive information of the mixing between the pn and the 
skin modes only from the isospin character of the tran- 
sition rates, although it is certainly useful for qualitative 
arguments. 

In this article we investigate character of the low- 
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energy El excitations via the RPA calculations on top 
of the Hartree-Fock (HF) wave functions. Calculations 
have been implemented for a number of doubly magic 
nuclei with various non-relativistic effective interactions, 
the Skyrme (SkM*, SkI2) [H 0, the Gogny (D1S, 
DIM) |28|,[29| and the M3Y-type (M3Y-P7) [Ia|3p[inter- 
actions, by using the numerical methods of Refs. [32h34| . 
Note that in these calculations the spurious cm. motion 
has been separated out to good precision, both for the 
Skyrme and the finite-range interactions. We then intro- 
duce a new method of decomposing the low-energy El 
strengths into the pn and the skin modes from the cal- 
culated transition densities, by which mixing of the two 
modes is taken into account in a straightforward man- 
ner. Via this decomposition we shall examine encrgy- 
and nucleus-dependence of character of the El transi- 
tions. Since interaction-dependence does not influence 
essential points of arguments, we shall mainly show re- 
sults of the SkI2 interaction below, with which observed 
PDRs in several nuclei are reproduced relatively well [35| . 

Transition density. The El transition operator is ex- 
pressed as 

(B1) = ^E^ (1) (^)-|E^ (1) ^): (1) 
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after the cm. correction. Here i is the index of nucleons 
and i £ p (i G n) indicates that the sum runs over protons 
(neutrons). The El strength from the ground state |0) 
to the 1~ state |a) in a even-even nucleus is B[E1; a) = 

2 

(a\0^ E1 ' |0) . As the lp-lh excitations, we calculate the 

transition amplitude {a\O iE1) |0) within the HF+RPA. 
By smearing out with Lorentzian, the strength function 
is given by 



7T ' 



1 



(u - uj a ) 2 + 7 2 
1 



(LU + LU a ) 2 + 7 2 



(a\O (E1) \0) . (2) 



We adopt 27 = 1 MeV throughout this paper. 

We next define the proton and neutron transition den- 
sity by 



it/) 



r 2 8pi x z =1) {r;a) = <J 



(for t z = p) 
(a\^2s(r-ri)riY^(ii)\0) 

(for r z = n) 



(3) 

Associated with in Eq. (TTJ, the El transition den- 

sity is defined as 

N 7 
6p( E1 \r;a) = -8p£= l \r-a) - -^ =1 )(r;a), (4) 
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which yields (a|C( B1 )|0) = / r 2 dr 5p {E1 ^ (r; a). It is no- 
ticed that, because the spurious cm. motion is well sep- 
arated, we have 



r 2 dr [<fp< A=1 > (r; a) + 5p£ =1) (r; a)] = . (5) 



The transition density Spi^ ^(r) is useful to investi- 
gate character of an individual transition. In the pn mode 
we have 5pp (r)+8pn (r) ~ at any r. In contrast, 
in the skin mode (5pp A_1 ' (r) and (r) have equal 

phases up to around the nuclear surface, and 5p^~^ (r) 
at larger r has the opposite phase to them, so that Eq. JSJ) 
should be fulfilled. This indicates that 5p^ E1 \r) is small 
in the interior because 5p^~ X \r) and 5p^n~ l \r) tend to 
cancel each other, while 8p^ E1 \r) (w (-Z/A) Sp { n =1) (r)) 
at large r mainly contributes to the El strength. 

In Fig. HJb) r 2 5 p^ -1 ^^] a) is depicted for several ex- 
cited states a in 132 Sn (also see Fig. [5] for the strength 
function). For reference, we show the proton and neu- 
tron density distributions in the HF ground state p Tz (r) 
in Fig. [Ha), along with the proton and neutron rms radii. 
The excitation to the uj a = 8.4 MeV state has typical 
transition densities of the skin mode, while the excitation 
to u> a = 12.7 MeV is dominated by the pn mode, although 
we may view skin-modc-likc behavior in r < 5 fm and in 

the displacement of the peaks between r 2 Spp X ~^ (r) and 
r 2 5pi^~^ (r). However, at energies in-between it is dif- 
ficult to classify the excitations into either of the pn or 
the skin mode. As the two modes seem to mix at these 
intermediate energies, the transition densities gradually 
change as oj increases. In this respect character change of 
the El strengths is not like a phase transition from the 
skin mode to the pn mode, should rather be regarded as 
crossover between the two modes. 

Decomposition of El strengths. It is reasonable to as- 
sume that El strengths are dominated by the two compo- 
nents, one corresponds to the skin mode while the other 
to the pn mode, up to around the GDR peak, though 
they can mix in individual excitations. In order to view 
the gradual change of the El character closely, we sep- 
arate the El transition density to each |a) into the pn- 
and the skin-mode components. Although it is not obvi- 
ous how the El strengths can be decomposed into these 
two ingredients, we here introduce a method which works 
practically well, particularly in handling their mixing. 
As argued above, while the pn-mode transitions hold IV 
character with 5pp ~~ *\r) and 8p^~ X ^ (r) having opposite 
phases, in the skin mode the transition densities have IS 
character in a certain spatial region, basically the interior 
of the nucleus, and are dominated by neutrons in the re- 
maining region. We therefore separate the excitation to 
each I a) into the pn- and the skin- mode components in 
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FIG. 1. (Color online) (a) Density distribution p p (r) (red 
line) and p n (r) (blue line) of the ground state in 132 Sn by the 
HF calculation with SkI2. The thick gray bar is bounded by 
the proton and neutron rms radii (r p and r n ), presenting their 
difference, (b) Transition densities r 2 <5pp A— ^ (r; a) (red line) 
and r 2 Spi^~^ (r; a) (blue line) to several 1~ states in 132 Sn 
by the HF+RPA calculation with SkI2. The overall phases 
are taken so as for the highest peak of r 2 5pn —1 \r; a) to be 
positive, for each a. 



a position-dependent manner by the following criterion: 



x ( Er , 

OPpn 

S E1 



(r;a) = 



for r satisfying 



Spy (r;a) 
bp n '(r;a) 



3pn ( 

*pS(r;a)=0 
otherwise . 



8p^ } (r;a)=5p( E1 \r;a) 



(6) 



Obviously Sp^ (r; a) = Sfffi (r; a) + 6p££ (r; a). To 
pick up components with (5pp A_1 '(r) w as the skin 
mode, it is appropriate to take < A s <C 1. For the 
12.7 MeV excitation of Fig. [TJ the transition density at 
r < 5fm and 9 < r < 13 fm is classified as Sp^iH ( r i a )- 
while the dominant part of 5 ^ r < 9 fm as 8p p ^ (r; a). 
Whereas it may be questioned whether the criterion of 



Eq. (|6]) is appropriate for very weak El excitations or 
El excitations at high energy, it will be useful in inves- 
tigating energy- dependent character of the El strengths 
below the GDR peak. Since the IS-like transition den- 
sity does not significantly contribute to the El strength, 
the condition of Eq. ^ indicates dominance of the neu- 
tron excitation for Sp^ irl (r), as long as transitions at 
relatively low energy are concerned. This separation of 
the spatial region according to the transition densities 
provides the way to decompose the El transition matrix 
element, 



(a\O^\0) = {a\O^\0) pn + {a\O^\0) skin : 
(a\O^\0) pn = j r^rSp^ina) 



{a\O^\0) BUin = / r*dr8p { £S(r-a) 



(7) 



and further to the decomposition of the strength func- 
tion, 

s^(«) = siT ) («) + sffl ) («) + s£ 1) («); 
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5 skin 
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(u + LU a ) 2 + 7 2 
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{a\O^\0) 



pn 



{lo -u a ) 2 +7 2 
1 



1 



(a|O (B1) |0) skin 



(lu - ui a ) 2 + 7 2 
1 



(u + LU a ) 2 + 7 2 



x2Rc 



t|0 (£1) |O}* (a|0( B1 >|O) skir 



The results of decomposing the El strength functions 
with X s = 0.05 are displayed in Fig. [2ta), along with 
B{E1) to an individual state \a), for 132 Sn. By this de- 
composition, we view more clearly the crossover from the 
skin mode to the pn mode in the El transition strengths. 
We have confirmed that the results do not change much 
with A s = 0.1 and A s = 0.2 [36[. Moreover, although 
the decomposition may reflect character of individual 
states | a), the decomposed strength functions do not 
have strong irregularity, after smoothing with respect to 
u>. 

We further define the ratios 



5 modc^i 



('mode' = 'pn', 'skin' or 'intf ') , 

(9) 

which are exemplified for 132 Sn in Fig. EJb). Note that 
-Rpn(w) = 1 — i?skin(^) — -Rintf(^)- These ratios are 
not sufficiently meaningful for ui at which S^ E1 >(u>) is 
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FIG. 2. (Color online) (a) Decomposed El strength func- 
tions and B(El;a) in 132 Sn. Solid line shows S {e1) (uj) ( = 



: (JS1) 



H + + s££V)), 



while blue-shaded and 

hatched areas present S^-^{u) and SSP(w), respectively. 
Red vertical bars display B(E1; a) without smoothing. Black 
arrow attached to the horizontal axis indicates the neutron 
threshold in the HF calculation, (b) R B ^ n (uj) and i?i n tf(w). 
Conventions correspond to those of the strength functions in 
(a). The SkI2 interaction is employed. 
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FIG. 3. (Color online) Decomposed El strength functions 
and B(El;a) with SkI2. See Fig. [2]for conventions. 



quite small, namely at very low u> (in absence of siz- 
able El strengths) as well as at high oj beyond GDR. 
However, above the energy at which i? s ki n (w) starts de- 
creasing (uj > 8 MeV in Fig. [5]), i? s kin(w) keeps decreas- 
ing almost monotonically up to the GDR peak energy, 
apart from small fluctuations. This illustrates that the 
position-dependent decomposition of the El strengths of 
Eq. ^ works practically well at u) < 20 MeV. 

The decomposed strength functions and the ratios are 
depicted for many doubly magic (or nearly doubly magic) 
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FIG. 4. (Color online) i? skin (u;) and i? int f(w) with SkI2. See 
Fig. [2] for conventions. 
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FIG. 5. (Color online) Comparison of ^^(cj) and _Ri nt f(ui) 
in 86 Ni among various effective interactions. See Fig. [2] for 
conventions. 



nuclei in Figs. [3] and |3J We confirm that the essential 
points found in 132 Sn hold for other nuclei. In neutron- 
rich nuclei there is a certain energy region where i? s ki n (w) 
exceeds 0.5, although i? s kin(w) values are not always quite 
meaningful at very low ui as in 132 Sn. Let us give special 
attention to the energy at which i? s ki n (w) comes down 
to 0.5 (after exceeding 0.5) in these neutron-rich nuclei. 
It is relatively high in the lighter nuclei 22 ' 24 0. How- 
ever, the energy at which i? s ki n (w) decreases down to 0.5 
is uj 10 MeV for all nuclei from 52 Ca to 132 Sn, with 
strikingly weak nucleus-dependence beyond the Ca re- 
gion. This consequence hardly depends on the effective 
interactions, since R s ^i n (uj) in the crossover energy re- 
gion is not quite sensitive to the effective interaction, as 
illustrated for 86 Ni in Fig. M 

PDR in 90 Zr. Because of the mixing of the pn mode 
and the skin mode, it is not guaranteed that all the ob- 
served low-energy El strengths should be attributed to 
the skin mode. Therefore careful study is needed in in- 
terpreting observed PDRs. We here argue character of 
the El strengths in 90 Zr as an example. 

The strength function and transition densities with 
SkI2 are shown in Fig. [6l in which a shoulder structure 
is found at w w 13 MeV. This shoulder state is a good 
pn mode but with two bumps in the transition densities. 
The other currently used interactions give similar results. 

Recent experiments have established a peak at uj ~ 
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FIG. 6. (Color online) (a) S {e1) {uj), B(El;a) and (b) 
r 2 6p (x=1) (r;a) in 90 Zr, with SkI2. See Figs. [1] and [2] for con- 
ventions. 



9 MeV in the El strength distribution of 90 Zr HQ. Sup- 
pose that the calculated strength at w w 13 MeV corre- 
sponds to this observed PDR, the PDR in 90 Zr is not 
associated with the skin mode, rather being a result of 
fragmentation of the pn mode. Despite the displacement 
in energy, it is still suggested that the observed PDR is 
not highly dominated by the skin mode, since mixing of 
the two modes is expected at u> ~ 9 MeV, according to 
nucleus-insensitive behavior of i? s km(w) in Fig. |4l Ex- 
perimental information on the transition density will be 
desired. 

Summary. By analyzing the transition densities ob- 
tained from the HF+RPA calculations in the doubly- 
magic nuclei, we have investigated character of the low- 
energy El excitations. In a certain energy region the 
pn mode (i.e. oscillation between protons and neutrons) 
and the skin mode (i.e. oscillation of the neutron skin 
against the inner core) mix with each other. We pro- 
pose a method to decompose the El excitations into the 
pn mode and the skin mode via the transition densities, 
which practically works well and enables us to handle 
their mixing explicitly. It has been shown that the El ex- 
citations gradually change their character as the energy 
increases. This behavior may be regarded as crossover 
transition from the skin mode at low energy to the pn 
mode at higher energy. The skin-mode strength comes 
down to a half of the full El strength at 10 sa 10 MeV from 
52 Ca to 132 Sn, generic for all nuclide and effective inter- 



actions under investigation as far as sizable El strengths 
are present at low energy 

Because of the energy-dependence of the character of 
the El excitations, it is not obvious whether observed 
PDR in an individual nucleus is dominated by the skin 
mode. We argue the case of 90 Zr as an example, via the 
calculated strength function as well as via the generic 
energy-dependence of the skin-mode fractions. We em- 
phasize importance of experimental information of the 
transition densities. 
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